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ABSTRACT 

We present a detailed analysis of the velocity distribution and orientation of orbits of 
subhaloes in high resolution cosmological simulations of dark matter haloes. We find 
a trend for substructure to preferentially revolve in the same direction as the sense of 
rotation of the host halo: there is an excess of prograde satellite haloes. Throughout 
our suite of nine host haloes (eight cluster sized objects and one galactic halo) there 
are on average 59% of the satellites corotating with the host. Even when including 
subhaloes out to five virial radii of the host, the signal still remains pointing out 
the relation of the signal to the infall pattern of subhaloes. However, the fraction of 
prograde satellites weakens to about 53% when observing the data along a (random) 
line-of-sight and deriving the distributions in a way an observer would infer them. This 
decrease in the observed prograde fraction has its origin in the technique used by the 
observer to determine the sense of rotation, which results in a possible misclassification 
of non-circular orbits. We conclude that the existence of subhaloes on corotating orbits 
is another prediction of the cold dark matter structure formation scenario, although 
there will be difficulties to verify it observationally. Since the galactic halo simulation 
gave the same result as the cluster-sized simulations, we assume that the fraction 
of prograde orbits is independent of the scale of the system, though more galactic 
simulations would be necessary to confirm this. 

Key words: methods: V-body simulations - methods: numerical - galaxies: forma¬ 
tion - galaxies: haloes 


1 INTRODUCTION 

There is mounting evidence that the Cold Dark Matter 
(CDM) structure formation scenario provides the most ac- 
cnrate description of our Universe. Observations point to¬ 
wards a “standard” ACDM universe comprised of 28% 
dark matter, 68% dark energy, and luminous baryonic mat- 
ter (i.e. galaxies, st ars, gas, and dust) at a mere 4% (cf. 
ISoergel et al. 1 1200,^1 . This so-called “concordance model” 
induces hierarchical structure formation whereby small ob¬ 
jects form first and subsequently merge to form progres¬ 
sively larger ob jects dWhite fc R.eesll978l : lDavis et al. Il985t 
iTormenlllQ^ . Whereas the large scale structure of our 
present Universe can be reconstructed very well by numer¬ 
ical simulations, the small scale structure still poses some 
problems. For example, there are many more subhaloes pre¬ 
dicted by cosmological s imulations than actually observed 
in ne arby galaxies (cf. iKlvnin et al. I IT999t iMoore et al. I 
^93)- The lack of observational evidence for these satel- 
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lites has led to the suggestion that they are completely (or 
almost completely) dark, with strongly suppressed star for¬ 
mation due to the removal of gas from the small proto¬ 
galaxies by the ionizing rad i ation from the first stars and 

S irs llBullock et alJDOnnl : iTullv et al. ll20o3 ISomerviIi3 
. Others suggest that perhaps low-mass satellites never 
formed in the predicted numbers in the first place, indi¬ 
cating problems with the ACDM model in general, replac¬ 
ing it with Warm Dark Matte r instead jCol t^tah 1 l2nonl: 
iBode. Ostriker fc Tiirokl 1200 ll : iKnebe et al. I 20oT R.ecent 
results from (strong) leasing statistics suggest that the pre¬ 
dicted excess of substructure is in fac t required to rec¬ 
oncil e some observation s with theory llPalal fc KochanekI 
I 2 OO 2 I: iDahle' et al. ]|20 o 3), although this conclusion has not 
been universally ac cepted llSchechter fc Wambsgansjl2002l : 

lEvans fc Wittl200.8l : ISand et al. ^2004^ . If. however, the leas¬ 

ing detection of halo substructure is correct and the over¬ 
abundant satellite population really does exist, it is vital to 
understand the orbital evolution of these objects and their 
deviation from the background dark matter distribution. In 
order to test the underlying ACDM model, more predictions 
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are necessary which can be confirmed or disproved by ob¬ 
servations. 

Here, we are investigating cosmological simulations 
based on the standard A-Cold Dark Matter model, con¬ 
cerning the sense of rotation of subhaloes. Host dark mat¬ 
ter haloes usually carry a small internal angular momen¬ 
tum, which is established by the transfer of ang ular mo¬ 
mentum from infalling matter via tidal t orques jPeeblesI 
1 19691 : iBarne^^^fetathiou ^98^. However, lGardnenl20(Hli 
as well as I'^witsk^^^ proposed another expla¬ 

nation for the origin of the angular momentum in galaxies 
and their dark matter haloes. They claim that haloes obtain 
their spin through the cumulative acquisition of angular mo¬ 
mentum from satellite accretion. These two descriptions are 
certainly linked together and mutually dependent, respec¬ 
tively. A detailed analysis of the orbits of satellite haloes 
shows that they are directly connected to the infall pattern 
of sa t ellites along the surrounding filaments (e.g. iTormenI 
I 1997 I : iKnebe et al. n i2004 IZentner et al. ]|20^. Those sub¬ 
haloes falling into the host at early times establish the angu- 
lar momentum of the inner regions of the primary halo (cf. 
IVitvitska et al. 1200211 and are channelled into the host along 
the same direction as those merging at later times. This leads 
to the speculation that satellites are preferentially corotat¬ 
ing with the host, which is the major assumption (and will 
be verified in) the current study. 

_ A similar study was recently presented bv lAzzaro et al. I 

i200,4l . They performed corresponding investigations us¬ 
ing both a cosmological simulation and observational data 
from the Sloan Digital Sky Survey. These observational data 
showed a signal of 61% corotating satellite galaxies, quite 
independent of magnitudes or distance to the host galaxy. 
When projecting their simulational data in order to “ob¬ 
serve” the simulations, they found a fraction of 55-60% pro¬ 
grade satellite galaxies. 

In another related study, IShaw et alJ (l2005^ studied a 
sample of 2200 (low resolution) dark matter haloes and de¬ 
termined the sense of rotation of all substructure particles 
with respect to the host halo. They found a very strong sig¬ 
nal for these particles to be corotating with the host. How¬ 
ever, we are using here a different method for determining 
the fraction of prograde satellites, classifying each satellite 
individually and then counting the number of pro- and ret¬ 
rograde satellites rather than (dark matter) particles. 

We present here the analysis of the sense of rotation 
for subhaloes in nine cosmological simulations: eight cluster¬ 
sized objects with varying merger histories and one galactic 
dark matter halo. These simulations are described in more 
detail in Section while Section |21 deals with the results 
when analyzing the full six dimensional phase-space infor¬ 
mation at hand. There we show how the sense of rotation 
for satellites can be determined and discuss various influ¬ 
ences on the fraction of prograde orbits. In Section 2] we 
investigate how our results can possibly be validated obser- 
vationally and conclude with a discussion and summary of 
our findings in Section Oandl^ 


2 THE SIMULATIONS 
2.1 The Raw Data 

Our analysis is based on a suite of nine high-resolution N- 
body simulations. Eight of them were carried out using the 
publicly available adaptive me sh refinement code mlapm 
jKnebe. Green fc Binnevll2nn focusing on the formation 
and evolution of a dark matter galaxy cluster containing 
of the order of one million particles, with mass resolution 
1.6 X 10® h~^ Mq and spatial force resolution ~2/i“^ kpc. 
We first created a set of four independent initial condi¬ 
tions at redshift 2 = 45 in a standard ACDM cosmology 
(Do = 0.3, Da = 0.7, Dt = 0.04, h = 0.7, ag = 0.9). 512® par¬ 
ticles were placed in a box of side length 64h“® Mpc giving a 
mass resolution of = 1.6 x 10® M©. For each of these 

initial conditions we iteratively collapsed eight adjacent par¬ 
ticles to one particle reducing our particle number to 128® 
particles. These lower mass resolution initial conditions were 
then evolved until 2 = 0. At 2 = 0, eight clusters (labeled 
C1-C8) from our simulation suite were selected, with masses 
in the r ange l-3x 10®'*h~^ M(T) and triaxiali ty parameters 
(cf. e.g. iFranx. Illingworth fc de Zeemj jl^^jUfrom 0.1-0.9 
(cf. Tablein. Then, as described by l'Hmnenrill997ll . for each 
cluster the particles within fiye times the yirial radius were 
tracked back to their Lagrangian positions at the initial red- 
shift (2 = 45). Those particles were then regenerated to their 
original mass resolution and positions, with the next layer 
of surrounding large particles regenerated only to one leyel 
(i.e. 8 times the original mass resolution), and the remaining 
particles were left 64 times more massiye than the particles 
resident with the host cluster. This conseryatiye criterion 
was selected in order to minimise contamination of the final 
high-resolution haloes with massiye particles. 

The ninth (re-)simulation was performed using the 
same (technical) appr oach but with the ART (Adaptiye Re - 
finement Tree) code iKraytsov. Klvpin fc Khokhlov^ll997^ . 
Moreoyer, this particular run (labeled Gl) describes the 
formation of a Milky Way type dark matter halo in a 
box of sidelength 20fe ~® Mpc. It is the s ame simulation as 
“Box20” presented in IPrada et al. I ^2^^5^ and for more de¬ 
tails we refer the reader to that study. The final object con¬ 
sists of about two million particles at a mass resolution of 
4 X Ilf h~^ Mq per particle and a spatial force resolution 
of 0.2h~^ kpc has been reached. All nine simulations haye 
the required resolution to accurately follow the formation 
and eyolution of subhaloes within their respective hosts and 
hence are well suited for the study presented here. 

The host haloes (as well as all substructure ob¬ 
jects down to 20 particles) are identified using the 
open source halo finder mhf^ (MLAPM’s-halo-finder; 
iGill. Knebe fc GibsorJl2004^ . mhf is based upon the adap¬ 
tive grid hierarchy of mlapm and acts with exactly the same 
resolution as the A-body code itself; it is therefore free of 
any bias and spurious mismatch between simulation data 
and halo finding precision arising from numerical effects. 
For every halo (either host or satellite) we calculate a suite 
of canonical properties based upon the particles within the 
virial/truncation radius. The virial radius Rvii is defined as 

^ MHF (and mlapm) can be downloaded from the following web 
page http://www.aip.de/People/aknebe/MLAPM 
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host 

Mvir 

Rvir 

A 

Vmax 

-^max 

*^tJ,host 

T 

age 

■^form 

/ AM \ 

\ AtM / 

^AMjM 

Cl 

2.9 

1355 

0.0157 

1141 

346 

1161 

0.365 

7.9 

1.052 

0.128 

0.125 

C2 

1.4 

1067 

0.0091 

909 

338 

933 

0.388 

6.9 

0.805 

0.122 

0.156 

C3 

1.1 

973 

0.0125 

828 

236 

831 

0.265 

6.9 

0.805 

0.100 

0.117 

C4 

1.4 

1061 

0.0402 

922 

165 

916 

0.639 

6.6 

0.750 

0.127 

0.207 

C5 

1.2 

1008 

0.0093 

841 

187 

848 

0.909 

6.0 

0.643 

0.129 

0.141 

C6 

1.4 

1065 

0.0359 

870 

216 

886 

0.073 

5.5 

0.567 

0.147 

0.153 

C7 

2.9 

1347 

0.0317 

1089 

508 

1182 

0.531 

4.6 

0.443 

0.844 

1.068 

C8 

3.1 

1379 

0.0231 

1053 

859 

1091 

0.587 

2.8 

0.237 

0.250 

0.225 

G1 

1.2x10-2 

214 

0.0229 

210 

44 

202 

0.491 

8.5 

1.232 

0.073 

0.107 


Table 1. Properties of the host haloes in our simulations. Masses are measured in Mq, velocities in km/s, distances in h~^ kpc, 

and the age in Gyr. (Mvir = virial mass, /?vir = virial radius, A = spin parameter, Vkiax = maximum of the rotation curve, /?max = 
position of the maximum, ciy host = velocity dispersion of the host, T = triaxiality parameter, age = time since half of the present day 
mass has formed, Zf^rm = corresponding formation redshift, ( aIai ) “ mean rate of relative mass change, — dispersion of the 

rate of relative mass change) 


the point where the density profile (measured in terms of 
the cosmological background density pb) drops below the 
virial overdensity Avir, i-e. M{< _Rvir)/(47r_R®iy3) = A virp6- 
This threshold Avir is based upon the dissipationless spher¬ 
ical top-hat collapse model and is a function of time for 
the given cosmological model. For 2 = 0, it amounts to 
Avir = 340. This prescription does no longer apply to sub¬ 
haloes where the point Rvii will not be reached due to the 
embedding of the satellite within the mass distribution of 
the host, i.e. the density profile will show a rise again at a 
certain point. In that case, we use this “upturn point” and 
truncate the object ignoring all particles outside of the cor¬ 
responding sphere. For a more elaborate discussion of this 
pr ocess and the halo finder, in particular, we refer the reader 
to lOill et al. I (|200^. 


2.2 The Host Haloes 


Our halo finder calculates a whole set of integral properties 
for each individual objec t, such as virial mas s Mvir, radius 
f?vir, spin parameter fcf. iBullock et ai~ll200lil 


A = 


L 

V^MviiVcRvlr: 


( 1 ) 


(with L = angular momentum, Vc = yTTMTir/Rfh = the 
circular velocity at Rvii, G = gravitational constant), the 
maximum of the rotation curve Vjnax (and its position i?max) 
and the eigenvalues a > fo > c of the inertia tensor, which in 
tur n can be used to construct the triaxia lity parameter (see 
e.g. iFranx. Illingworth fc de Zeeuw^ll991^ 


T = 



( 2 ) 


In order to get a quantitative measure for the mass ac¬ 
cretion history for each of our host haloes, we also compute 
the dispersion of the fractional mass change rate; 


_2 _ 1 f AMi / AM 

^AM/M - ^ 2^ [AtiMi ~\AtM/) 

i=l 


( 3 ) 


where Nout is the number of available outputs from forma¬ 
tion 2form to redshift 2 = 0, AM; = |M( 2 i) — M( 2 i+i)| the 


change in the mass of the host halo, and AU is the respec¬ 
tive change in time. The formation redshift 2form is defined 
as the redshift where the h alo contains half of its present day 
mass llLacev fc Colell993^ and determines the age of the ob¬ 
ject. The mean growth rate for a given halo is calculated as 
follows 

/^\ _ 1 _ AMi , , 

\ AtM / Aout - 1 ^ AtiMi ■ ^ 

i=l 

A large dispersion (Jam/m now indicates a violent formation 
history whereas low values correspond to a more quiescent 
formation. As our definition for formation time implies that 
the host halo’s mass grows by a factor of two until 2 = 0, 
an estimate for the mean growth rate is simply given by the 
inverse of the host’s age and should be compared against 
the rate calculated via Eq. 0. We summarize these values 
alongside other relevant properties in Tabled 


3 ANALYZING THE SIMULATIONS 

3.1 Defining the sense of rotation 

In order to determine whether a satellite rotates in the same 
direction as the host halo (prograde orbit) or in the opposite 
direction (retrograde orbit), we need to evaluate the align¬ 
ment of the host’s rotation axis and the orbital rotation axis 
of the satellites. 

The rotation axis of objects, such as host dark matter 
haloes, is not well defined^, thus we use the angular momen¬ 
tum for the description of its rotation: 

L{R) = ^ mifi X Vi , (5) 

ri<R 

which is being calculated at various distances R from the 
host’s centre using all interior particles. The (dark matter 
particle) velocities Vi here are in the rest frame of the host 
as are the particles’ positions f). 


^ For a rigid body, the relation between the angular velocity vec¬ 
tor uj determining the rotation axis and the angular momentum 
vector L can be derived via L = lul, where I is the inertia tensor. 
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Figure 1. The distribution of directions of orbital angular mo¬ 
menta on a unit sphere. The host’s internal angular momentum 
points from the centre to the North Pole. The points where the or¬ 
bital angular momenta of satellites pierce through the unit sphere 
are marked with a More points in the Northern hemisphere 
than in the Southern hemisphere indicate a majority of prograde 
orbits. 


for C4 clearly shows a majority of points in the Northern 
hemisphere, indicating an excess of prograde orbits. In con¬ 
trast, the distributions for the hosts C3, C5 and C8 appear 
highly isotropic. A quantitative analysis of these hrst im¬ 
pressions will follow in the next sections. 


3.2.1 Differential distribution of angles 

We determined the cosine of the angle between the inner 
host halo angular momentum L 30 and the orbital angu¬ 
lar momentum L°at for every satellite residing within the 
virial radius i?vir of the host halo. Figure |21 now shows the 
(binned) differential distribution of cos(j). We readily ob¬ 
serve a (slight) increase in the fraction of satellites for larger 
values of cos (j>, indicating an overbalance of prograde orbits. 
While the solid lines represent the distributions when count¬ 
ing each satellite once (“number weighing”), the dotted line 
is derived by weighing the contribution of each satellite by its 
mass. Comparing both curves demonstrates the importance 
of high-mass satellites on the prograde signal. Unfortunately 
the trend is not that obvious and we rather conclude that 
the fraction of prograde subhaloes is not significantly influ¬ 
enced by heavier satellites. The red solid horizontal line in 
each plot marks the expected distribution for an isotropic 
distribution of cos (f. 


In order to eliminate the influence of (massive) sub¬ 
structures at large halocentric distances we used the an¬ 
gular momentum vector L30 as defined by the material 
inside a sphere containing 30% of the host’s virial mass 
(« 0.2i?vir on average). This also ensures to focus on the 
pro perties of the inner halo whose material collapsed hrst 
fcf. lHelmi. White fc SDril^l2^03^ . 

The rotation axis of a satellite orbit simply coincides 
with its orbital angular momentum vector 

Lsat ~ nisatrsat X Vsat ■ (6) 


where fsat and Vsat are again measured in the rest frame of 
the host halo. 

To check whether a satellite is co- or counterrotating we 
use the angle between the internal angular momentum of the 
host and the orbital angular momentum of the satellite, i.e. 
more precisely we use the scalar vector product 


cos (f 


f* ?orb 
-^30 • -^sat 


(7) 


and dehne a satellite to be prograde for cos (j> > 0 and retro¬ 
grade for cos (f> < 0. Thus, we also include satellites on nearly 
“perpendicular” orbits (cos (j> ~ 0), whose classihcation as ei¬ 
ther pro- or retrograde may be questionable. However, such 
subhaloes are part of our sample and any hypothetical ob¬ 
server of the system will have to deal with them, too. We 
will come back to this point in Section 2] where we project 
our data into the “observer’s plane”. 


3.2 Determining the sense of rotation 

For a visual impression of the distributions of orientations 
of subhaloes and their (an-)isotropy, we plot the direction 
of the orbital angular momenta on a unit sphere in Figure^ 
(projected to a plane using the Aitoff-Projection). The panel 


3.2.2 Cumulative distribution 

Presenting a differential distribution always includes the 
problem of binning the data and hence introducing (unnec¬ 
essary) noise, especially for small number statistics. To avoid 
this and provide a smoother view of the data, we examine the 
cumulative distribution function, as shown in Figure 0 Here 
we present the number fraction (and the mass fraction, dot¬ 
ted line) of satellites with cos (p smaller than a given value. 
For an isotropic distribution we expect a straight diagonal 
line as indicated by the red solid line. We can see, that the 
distributions in Figure|21do not seem to be isotropic, just as 
already indicated by Figure]^ However, this does not nec¬ 
essarily lead to a bias towards prograde orbits. If we still 
had the same number of pro- and retrograde orbits, the dis¬ 
tributions would be point symmetric. But we observe the 
clear trend in (most of) our systems that the data curves lie 
below the isotropic expectation. 


3.2.3 Confirming the signal via a Kolmogorov-Smirnov 
test 

In order to validate the preference for satellites to be on 
prograde orbits as seen in both Figsandand to quantify 
the signihcance of the si gnal, we carried ou t a Kolmogorov- 
Smirnov (KS) test le.g. IPress et al. ]|20 o3 ) for the number 
weighted distributions. The KS test provides us with the 
maximum deviation of the data set from the isotropic dis¬ 
tribution {D) and the probability (p) for the data being 
consistent with isotropy. These numbers are listed in Ta¬ 
ble 1^ and show that only the distributions of orientations 
in simulations C3 and C5 (and possibly C8) are close to 
isotropy with a nearly identical number of pro- (rip) and 
retrograde (wr) orbits. The latter numbers are also given in 
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Figure 2. The distribution of the cosines of angles cos </< between the orbital angular momentum of the satellites and the internal angular 
momentum of the host halo for nine different simulations (first panel = oldest halo) at redshift z = 0.0. Shown is the fraction of the 
number of satellites per cos(/)-bin (black solid line), divided by the size of each bin. The red horizontal line marks the value for isotropic 
distribution. The dotted curve represents the corresponding mass fraction of the satellites per bin. 


Table 121 supplemented by the fraction of subhaloes on pro¬ 
grade orbits, i.e. fp = npjijip -frir). In all other simulations 
we find a quantitative conhrmation of an excess of prograde 
orbits. 

However, as mentioned above, a deviation from isotropy 
does not automatically mean an overbalance of prograde or¬ 
bits: it is the asymmetry of the distribution, which could 
validate the result. If the satellites were equally distributed 
between pro- and retrograde orbits, the differential distribu¬ 
tion would be symmetric about the vertical line cosi^ = 0 
and the cumulative distribution function would be point 
symmetric around (0, /(< 0.5)). We can test for this (point) 
symmetry by plotting the cumulative distribution not in the 
interval [—1,-|-1] but rather in reverse with cos<)) ranging 
from [-1-1,—1]. If these two ways of calculating the distri¬ 
bution give rise to identical curves, then the distribution is 
(point) symmetric, without any preference for pro- or retro¬ 
grade orbits. In order to further quantify this kind of sym¬ 
metry, we can again apply a KS test to these two “mirror”- 
distributions. For six of our nine simulations we could not 
detect any point symmetry, the probability for point symme¬ 
try is at most 5%. Only for the simulations C3 (84.4%), C5 
(88.5%) and C8 (19.8%) we do find a stronger sign of sym¬ 


metry. These systems are the cases, which also show more 
or less isotropic distributions and thus this result comes at 
no surprise as isotropy entails (point) symmetry. We will 
elaborate upon these three systems later on in Section roi 
However, for all the other simulations, this analysis confirms 
the credibility of our results which can be summarized by the 
mean fraction of prograde orbits of 59% ± 7% throughout 
all simulations. 

The strength of the signal can be measured not only 
by the (number) fraction of prograde orbits as given in the 
last column of Table |21 but also by (the cosine of) the angle 
between the host’s angular momentum and the net angular 
momentum of all subhaloes 


cos $net 


7 V^^sat r*orb 
-^30 • -^sat,i 

Laol Erir kti\' 


( 8 ) 


These values are also listed in Table|21and with one exception 
(C5, the most isotropic case) we only have positive values, 
stressing again the preference of prograde orbits. One needs 
to bear in mind though that the use of cos 4>net rather con¬ 
firms a “mass weighted signal” as the angular momentum 
^saM depends on the mass of each individual subhalo. 
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COS0 COS0 COS0 


Figure 3. The cumulative distribution of the cosines of angles cos <f> between the orbital angular momentum of the satellites and the 
spin of the host halo for nine different simulations (first panel = oldest halo) at redshift z=0.0 (black solid line). The dotted line shows 
the cumulative distribution with mass weighing - instead of just counting the satellites, the masses of the satellites were added. A strong 
increase in that curve indicates a massive satellite. The red line shows the ideal case of an isotropic distribution. 


3 . 2.4 Uncertainties in the host’s rotation axis 

For the analysis so far, we used the angular momentum of 
the inner halo, i.e. L30 as defined by the material inside a 
sphere encompassing 30% of the virial mass. We now ex¬ 
amine the effect when using the angular momentum based 
upon larger fractions, i.e. 50% or even 100%, of the host’s 
virial mass. 

We expect variations in the host’s angular momentum 
for two reasons: first, a dark matter halo is not a rigid body 
and secondly, the angular momentum in the outer parts may 
be “contaminated” by recently accreted material and in¬ 
falling satellites. These variations in the direction of the host 
angular momentum vector compared to the angular momen¬ 
tum of the inner 30% are quantified in Table |21 where we list 
the angle between the two respective vectors 


0 , 50/100 


= COS ^ 


( L30 ■ T50/100 
^ 30 ^^ 50/100 


(9) 


The last line in Table summarizes the mean angle when 
averaging over all hosts along with its standard deviation. 

More interesting now are the variations in the fraction 
of prograde orbits, i.e. A/p = /p — /p°. Table 31 summa¬ 


rizes the corresponding prograde fractions (/p° and /p°°) 
and their deviations A/p = fp — fp^ from our reference “in¬ 
ner” fractions /p°. Even though there are minor variations 
amongst individual haloes, the mean fraction of prograde 
orbits stays rather constant irrespective of the definition for 
the host’s angular momentum. 


This result can be understood considering the fact that 
there are not many satellites on perpendicular orbits as in¬ 
dicated by Figure Q- the orbital angular momenta prefer¬ 
entially point to the region around the poles rather than to 
the equator region. Since we define all satellites as corotat¬ 
ing whose angular momenta lie north of the equator, varia¬ 
tions in the direction of the angular momenta do not have 
a prominent influence on the number of prograde orbits. 


Hence, despite the rather large differences in the relative 
orientation of the respective angular momenta (cf. Tabled 
the influence of the point, where to define the host’s angular 
momentum, on the mean fraction of prograde orbits is only 
minor. 
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host 

D 

P 

cos illnet 

^^sat 

Up 

Ur 

fp 

Cl 

0.120 

0.064 

0.959 

233 

134 

99 

0.575 

C2 

0.113 

0.387 

0.318 

124 

72 

52 

0.581 

C3 

0.073 

0.991 

0.248 

68 

37 

31 

0.544 

C4 

0.220 

0.008 

0.897 

109 

77 

32 

0.706 

C5 

0.062 

0.990 

-0.286 

97 

46 

51 

0.474 

C6 

0.198 

0.027 

0.978 

106 

73 

33 

0.689 

C7 

0.124 

0.050 

0.369 

234 

145 

89 

0.620 

C8 

0.059 

0.617 

0.206 

322 

171 

151 

0.531 

G1 

0.156 

10-® 

0.622 

557 

331 

226 

0.594 





< fp 

> = 

0.590 ± 0.070 


Table 2. Isotropy and prograde fraction. The isotropy was 
checked with a KS test, which calculates the maximum devia¬ 
tion D and the probability p for the consistency with an isotropic 
distribution, cos $net is the cosine of the angle between the in¬ 
ner host angular momentum and the net angular momentum of 
all satellite orbits, rigat is the total number of satellites that can 


host 

f50 

Jp 

(A/f) 

floo 

Jp 

(A/i°0) 

Cl 

0.601 

(±0.026) 

0.579 

(±0.004) 

C2 

0.556 

(-0.024) 

0.605 

(±0.024) 

C3 

0.471 

(-0.074) 

0.515 

(-0.029) 

C4 

0.706 

(±0.000) 

0.725 

(±0.018) 

C5 

0.588 

(±0.113) 

0.660 

(±0.186) 

C6 

0.689 

(±0.000) 

0.585 

(-0.103) 

C7 

0.585 

(-0.034) 

0.546 

(-0.073) 

C8 

0.481 

(-0.050) 

0.488 

(-0.043) 

Cl 

0.639 

(±0.045) 

0.627 

(±0.032) 

< fp > 

0.591 ±0.077 

0.592 ± 0.069 


Table 4. Quantifying the influence of the host’s angular mo¬ 
mentum (at 50% and 100% of the host mass) on the fraction of 
prograde orbits. The last line gives the mean prograde fractions 
with standard deviations. 


be split into rip and rir, the number of prograde and retrograde 
orbits, respectively. The prograde fraction fp is given in the last 
column and the last line simply gives the mean fraction (along 
with its standard deviation) of prograde orbits when averaging 
over all hosts. 


host [°] alOO [°] 


Cl 

32.6 

27.0 

C2 

19.1 

46.5 

C3 

28.4 

35.4 

C4 

2.4 

7.4 

C5 

48.8 

37.3 

C6 

3.4 

17.4 

C7 

40.5 

43.0 

C8 

16.9 

65.2 

Cl 

15.3 

32.8 


<o> 23 ± 15 35 ±17 


host 

1.5 

riscit 

-Rvir 

fp 

3.0 ■ 

risat 

Rvir 

fp 

5.0 ■ 

nsat 

Rvir 

fp 

Cl 

458 

0.544 

750 

0.507 

774 

0.509 

C2 

243 

0.560 

384 

0.534 

411 

0.530 

C3 

133 

0.534 

192 

0.599 

201 

0.602 

C4 

212 

0.693 

311 

0.695 

318 

0.686 

C5 

144 

0.472 

210 

0.471 

226 

0.469 

C6 

199 

0.638 

365 

0.592 

395 

0.587 

C7 

382 

0.613 

632 

0.568 

659 

0.561 

C8 

446 

0.545 

570 

0.549 

604 

0.558 

Cl 

904 

0.579 

1589 

0.566 

3015 

0.512 

< fp > 

0.575 ± 0.061 

0.565 ±0.060 

0.557 ± 0.060 


Table 5. Numbers of satellites and prograde fractions when in¬ 
cluding satellites from the immediate environment, i.e. including 
satellites outside the virial radius of the host halo, up to 1.5, 3.0 
and 5.0i?vir- The mean values of the prograde fractions < fp > 
for all simulations are provided in the last line. 


Table 3. The difference in the direction of the host angular mo¬ 
mentum at radii including 50% and 100% of the host mass, with 
respect to the 30%-angular momentum. 


stantially, the sense of orientation of the satellites is hardly 
affected. The signal weakens becoming more isotropic, but 
yet remains even at a distance of bRvir- 


3.2.5 Environmental dependence 

The recent study of IPrada et al. I j2Q05h showed that iso¬ 
lated galaxy-sized haloes display all the properties of relaxed 
objects up to 2-3i7vir- One of their haloes under investiga¬ 
tion (i.e. “Box20”) w as in fact the galactic halo a lso used 
in our study. Further, Em Knehe fc GibsorJ ^200!^^ showed 
that there exists a prominent population of backsplash sub¬ 
haloes: satellites that once passed the host’s virial radius 
but are now residing in the outskirts of the halo. It there¬ 
fore appears interesting to investigate the effect of allowing 
for satellites outside the virial radius to be included in the 
determination of the prograde fraction. The results for grad¬ 
ually increasing the trading area and considering satellites 
within a sphere of radius 1.5, 3, and 5i?vir are given in Ta¬ 
ble 1^1 Besides listing the updated prograde fraction we also 
present the number of satellites in the respective sphere. 
We can see that while the number of subhaloes grows sub- 


3.3 Discussion 

We have seen in our simulations that the average fraction of 
prograde orbits is 0.590 ± 0.070. 

In Section 0 we sketched a scenario, where matter 
streams into the (progenitors of) host haloes from preferred 
directions, i.e. along the surrounding filaments, thus estab¬ 
lishing the angular momentum of the host. Since today’s 
satellites are also expected to fall into the host from the same 
preferred directions, this picture suggests an overabundance 
of prograde orbits. This scenario is supported by the anal¬ 
ysis presented in the previous Section Even satellites 

as far away from the host as 5 Rvii will eventually fall into 
the host in a way that complies with the already established 
sense of rotation. 

Though this picture in general seems to be confirmed, 
we still would like to understand why not all our systems do 
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exhibit the signal. If we use the information about the hosts 
presented in Tableland compare their integral properties 
with the prograde fraction for each halo, we note a slight 
trend of larger prograde fractions with increasing spin pa¬ 
rameter (neglecting the galactic halo for the moment). The 
isotropic hosts C3 and C5 possess only a very small spin 
parameter of A = 0.0125 and 0.0093, respectively. This con¬ 
nection is understandable when considering the origin of the 
host angnlar momentnm from mergers with subhaloes - the 
smaller the angular momentum, the more likely its direction 
is “switched” after yet another merger. 

However, host halo C8 seems to be the exception to this 
rule: just like C3 and C5 it is one of the few hosts with a very 
modest tendency for corotating satellites, but has an average 
spin parameter (A = 0.0231). However we need to keep in 
mind that it is the youngest system under investigation. And 
its evolution gives us the key to understanding the missing 
signal: from a visual inspection of its mass accretion history 
and particle distribution over time we could infer that it 
just recently merged with two other halo es of nearly the 
same size, i.e. it is a (recent) tr iple merger! Ivitvitska et al. I 
^2002^ (see also lGardne3l200lh found in their simulations 
that haloes with recent mergers (after z = Z) generally show 
a larger spin parameter than “quiet” haloes. But even if the 
internal angular momentum of the resulting host halo has 
a “stable” direction, the satellites have been swirled around 
violently during the merging event and thus can be expected 
to have rather random orbit orientations, just as we found 
in the case of C8. 


4 “OBSERVING” THE SIMULATIONS 

We have seen in our simulations, that subhaloes are pref¬ 
erentially moving on prograde orbits - around 60% of the 
subhaloes are corotating with the host halo. Could this sig¬ 
nal be observed? Assuming that all our satellite haloes carry 
enough baryonic matter, in particular stars and gas, to be 
detectable by an observer, we are facing the question how 
the systems in our simulations appear from the observer’s 
point of view. 

4.1 Cumulative line-of-sight velocity distribution 

Usually an observer is not in the fortunate situation to know 
the positions and velocities of the satellites in all three di¬ 
mensions. Yet it is possible to determine the velocity of a 
satellite along the line of sight and its projected position rel¬ 
ative to the (density) centre of the host halo. In order to get 
a first impression of the velocity distribution of satellites, we 
looked at the cumulative distribution of line-of-sight veloc¬ 
ities with respect to the host halo. We chose one hundred 
(isotropically distributed) random lines of sight and then 
averaged over the individual ID velocity distributions. The 
resulting mean line-of-sight velocity distributions (with re¬ 
spect to the rest frame of the host halo and normalized to 
the velocity dispersion cr„, host) are shown in Figured The 
error bars represent the standard deviation. 

If the velocities of the satellites were isotropically dis¬ 
tributed (i.e. random directions with respect to the line- 
of-sight, but the same absolute values), we would observe 
the thick red lines in Figure 2] Our data curves (thin black 



Figure 5. Classifying observed satellites into prograde and retro¬ 
grade orbits: The figure shows the projection plane perpendicular 
to the line of sight. Lp is the projected angular momentum of the 
host halo, Arp is the projected distance between the host halo 
centre and a satellite (marked with an asterisk). The angle a be¬ 
tween this distance vector and the vector from the host halo centre 
to the red end of the host (fy^u) is used for deciding whether the 
satellite is located on the red end or the blue end side of the host 
halo. 

lines) are coinciding nearly perfectly with their correspond¬ 
ing isotropic curves. A KS test simply confirms the eye- 
balled result: the probability of the data curves being con¬ 
sistent with the isotropic curves is practically 100%. This 
is actually not too surprising as the signal for a majority of 
prograde orbiting satellites is rather weak, but we could have 
expected a scattering in the data curves introduced by ob¬ 
serving the halo from various viewing angles. However, the 
standard deviations for data and isotropic curves are nearly 
identical and rather small. The line-of-sight velocities are 
isotropically distributed, and hence we need to devise a dif¬ 
ferent approach to observationally confirm the existence of 
prograde orbits. 

4.2 Classifying prograde and retrograde orbits 

In what follows we lay out a method allowing an observer to 
distinguish prograde and retrograde orbiting satellites (cf. 
Figure |3 which has already been used bW^zar^^^l/ 

I 2005) in a similar analysis (see also e.g. IZaritskv et al. 

IQQTT , for an observational study): One first decides 
whether the satellite is moving away from or towards the 
host by determining the line-of-sight velocity Auios with re¬ 
spect to the host halo. This needs to be compared with the 
host’s sense of rotation: We call the approaching side of the 
host the “blue end”, while the receding part is called “red 
end” and define 5=1 and S' = — 1, respectively. Based upon 
this information we can now define a (one-dimensional) pe¬ 
culiar velocity for each satellite 

Upec — S • Auios , (^^) 

and hence a satellite is on a prograde orbit for Upec> 0 and 
retrograde for Upec< 0 . 

4.3 Peculiar velocity distributions 

In Figure 1^ we present the results for the (differential) pe¬ 
culiar velocity distribution. We used again one hundred ran¬ 
dom lines of sight and the velocities are normalized by the 
hosts’ velocity dispersion. All satellites to the left of the 
vertical dash-dotted “zero” line have been classified as ret¬ 
rograde, while all satellites to the right are observed to be 
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Figure 4. The cumulative distribution of line-of-sight velocities, normalized to the hosts’ velocity dispersion (t„ host a-nd averaged over 
the results from one hundred random lines of sight. The thick red lines are ideal distributions for isotropic directions of velocities (but 
the same absolute values). Error bars represent the standard deviation for a few sample points. There is practically no difference between 
the isotropic and “observed” curves, indicating an isotropic distribution of line-of-sight velocities. 


prograde. In addition, the mean number of prograde and 
retrograde orbits (along with the la deviation) is given for 
every host halo in the upper right corner of every panel. The 
signal of an excess of prograde orbits as found when using the 
full three-dimensional (velocity) information is confirmed, 
but the signal is weakened in comparison to our previous 
findings. We only get a mean observed prograde fraction of 
53.3% (±3.3%) as opposed to 59.0% from Section !!). 2. 31 The 
individual observed prograde fractions for each dark 
matter host are given in the inset panels of the respective 
plot of Figure!^ 


4.4 Discussion of the “observer’s” result 
4-4-1 Misclassification of orbits 

The weakening of the prograde fraction for the observer 
is the result of a possible false classification of the satel¬ 
lites’ sense of rotation. The method applied (i.e. Eq. innt ) 
is based upon one not yet mentioned assumption: it will 
only work faultlessly, if the satellites are on circular orbits, 
which usually is not the case. Previous studies have shown 
that satellites can be on highly eccentric orbits, and their 


eccentricity distributions can be fit ted by a Gaussian with 
a mean eccentric i ty of ep = 0.61 _jGilFe^lJ_j2004[j_^t 
see also iTorm enI Jl997l) ; lTormen'~ni^fen^^^TOTl Il99fji : 
iGhigna et al. I lll99Sh ). 

In observational studies (e.g. Zaritskv et al. I ll 9971 : 
ICarignan et al. Ill997t lAzzaro et al. I 2005ir ~this method is 
used without mentioning its limitations, thus we here ex¬ 
plain the cause of misclassification in more detail. 

Figure Q elucidates the situation by viewing a sample 
satellite-host system from “above”, i.e. perpendicular to the 
line-of-sight, which lies in the plane of the paper. In compli¬ 
ance with our terminology, the part of the host halo which 
is hashed with red lines is called “red end” - it is the part 
of the host halo that is receding from us. The “blue end” is 
defined accordingly. The ellipse represents a (idealised) part 
of an eccentric satellite orbit whose measured relative line- 
of-sight velocity Auios is positive in the red dotted part and 
negative in the blue dashed part. The small velocity vectors 
attached to the ^-symbol for the satellite at the points 1, 
2 and 3 should indicate that the satellite is corotating with 
the host halo and hence on a prograde orbit. However, we 
can readily see in Figure 0 that the satellite will be mis- 
classified, if it happens to be observed in the blue (dashed) 
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Figure 6. The peculiar velocity distribution (normalized to the hosts’ velocity dispersion host) for nine simulations. Here, velocities 
are defined to be positive, if the observer classifies the satellite as prograde, negative velocities mean retrograde motion. 


region of its orbit but on the red end side of the host halo 
(cf. positions 2 and 3). Even though the satellite shares the 
sense of rotation with the host, it appears to be counter¬ 
rotating; it approaches us while the host is receding. Such 
situations can occur when the turning points of the orbit, 
i.e. the points where the line-of-sight velocity of the satellite 
changes its sign [F is the farthest and N is the nearest point 
to the observer), are not aligned with the line of sight. 


4-4-2 Origin of misclassification 

In order to support our theory of misclassification, we con¬ 
sider its sources in more detail. We are in the unique sit¬ 
uation where we can distinguish between “correctly” and 
“falsely classified” orbits as we can compare the observer’s 
classification to the three-dimensional classification used in 
Section 0 

We argued that the non-circularity of the orbits is the 
(possible) origin for a misclassification. In order to verify this 
conjecture, we use the circularity parameter Acirc, defined as 
the ratio of the angular momentum of the satellite orbit Lsat 
to the angular momentum of a correspon ding circular orbi t 
Lcirc with the same orbital energy Esat icf. lGill et al. l2004l : 


1 ffsat (f-^'sat) /1 1 \ 

= J 7 (ft) 

i^circ\-C/sat ) 

Figure |H1 shows the binned circularity distribution for all 
satellites (black line), correctly classified satellites (blue 
hashed area) and falsely classified satellites (red hashed 
area). Satellites with circularity Acirc = 1 are on circular 
orbits whereas radial orbits will have Acirc = 0. However, 
a circularity of, for instance, A circ = 0-9 corresp onds to an 
eccentricity of 0.6 (cf. Fig.11 in IfliH et al. 11200^ and hence 
small deviations from Acirc = 1 already imply rather ec¬ 
centric orbits! Figure (Whence proves that satellites with a 
high circularity are more likely to be correctly classified than 
their corresponding low circularity counterparts, stressing 
that the deviation from circular orbits increases the proba¬ 
bility for misclassification. 


4-4-3 Quantitatively understanding misclassifications 

In Appendix we show how the observed prograde frac¬ 
tion will be weakened due to misclassification, assuming that 
there is an equal probability of wrong classification for both 
pro- and retrograde orbits (the variation in our simulations 
is at most a few percent). This follows from the simple con¬ 
sideration that all wrongly classified prograde orbits will ac- 
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host fp c /*= 


Cl 

0.575 

0.665 

± 

0.060 

0.525 

± 

0.009 

0.526 

± 

0.025 

C2 

0.581 

0.656 

± 

0.077 

0.525 

± 

0.012 

0.525 

± 

0.032 

C3 

0.544 

0.665 

± 

0.091 

0.515 

± 

0.008 

0.506 

± 

0.041 

C4 

0.706 

0.684 

± 

0.067 

0.576 

± 

0.028 

0.609 

± 

0.027 

C5 

0.474 

0.666 

± 

0.081 

0.491 

± 

0.004 

0.497 

± 

0.035 

C6 

0.689 

0.638 

± 

0.075 

0.552 

± 

0.028 

0.546 

± 

0.036 

C7 

0.620 

0.679 

± 

0.082 

0.543 

± 

0.020 

0.530 

± 

0.051 

C8 

0.531 

0.676 

± 

0.060 

0.511 

± 

0.004 

0.499 

± 

0.014 

G1 

0.594 

0.699 

± 

0.058 

0.538 

± 

0.011 

0.555 

± 

0.023 

< • > 

0.590 

0.670 

± 

0.016 

0.531 

± 

0.024 

0.533 

± 

0.033 


Table 6. A summary of prograde fractions: fp is the fraction derived using the 3D information provided by the simulations (cf. Tahlel^. 

the expected observed fraction of prograde satellites based upon Eq. ini and fp^^ the actual observed fraction (cf. Figure IHl. c is 
the mean fraction of correctly classified orbits (from one hundred lines of sight). Mean values, averaged over all simulations, are given 
again in the last line. 



line of sight 


Figure 7. The host halo and one of its satellites on an idealised 
elliptic orbit as seen from “above”. The observer is somewhere 
beyond the bottom of the page, looking at this system along the 
line of sight. HC marks the host halo centre, F is the farthest 
and N is the nearest point of the orbit relative to the observer. 
The part of the orbit when the satellite is approaching is dashed 
(blue), the receding part is dotted (red). Intersections of the orbit 
with the imaginary line dividing the host halo into red and blue 
end are named Sp and Sf<[. The satellite actually is on a prograde 
orbit, but will be misclassified as retrograde, if it is observed at, 
for instance, position 2 or 3 (between F and Sp or Spi and A”, 
respectively.) 



0.0 0,2 0,4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
circulorlty clrculorlty clrculority 


Figure 8. Circularity distribution of satellites, i.e. the fraction 
of satellites as a function of circularity Adrc divided by the size of 
each bin (solid lines). For circular orbits we have a higher fraction 
of correctly classified satellites (blue hashed area) than misclassi¬ 
fied satellites (red hashed area). 


tually be categorized as retrograde and vice versa. If there 
is, for instance, a majority of prograde satellites, more satel¬ 
lites will “move” from prograde to retrograde than from 
retrograde to prograde and hence the prograde fraction is 
weakened. We yield the following equation for the observed 
ratio of prograde Up to retrograde orbits 


obs 

g 


^ _ 

oobs 


c - q + [1- c) 
c-\- q- (1 — c) 


( 12 ) 


with q = Up/Ur being the ratio of prograde to retrograde 
orbits in the 3D-simulation and 


7~ip TLj' 


(13) 


being the probability of right classification. 

A majority of prograde orbits in our simulations can 
only change into a minority for an observer, if c becomes 
less than 50%, provided that our statistical sample is big 
enough. We determined the mean fraction of correct clas- 
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Figure 9. The fraction of correctly classified satellites as a func¬ 
tion of relative line-of-sight velocity Auios- The probability for 
correct classification is higher for high velocities. The error bars 
are the standard deviation from averaging the results over one 
hundred random lines of sight. 

sification c for every simulation by averaging over hundred 
randomly chosen lines of sight. We found c to be around 
67%, the individual values are listed in Table [SI From this 
value for c and the ratio q taken from our 3D-simulations, 
we can now use Eq. iim to deduce an expected observed 
prograde fraction These predicted prograde fractions 

are given in Table ini and can be compared to the actually 
observed fractions (cf. Figure]^. We hnd that our antic¬ 
ipated observed prograde fractions fit very well the mean ob¬ 
served prograde fractions. The mean value for the observed 
prograde fraction of all simulations was found to be 53.3%, 
which is close to the predicted value of 53.1%. 

4-4-4 Overcoming misclassification 

The observer usually does not know the circularity of the or¬ 
bit and also has no other means to check whether the satel¬ 
lite is residing in a disadvantageous region or not. However, 
he will be able to determine, if the satellite is very close to 
the border of such a region: close to the turning points (the 
nearest or farthest orbit point, V or F in Figure |7| respec¬ 
tively) its line-of-sight velocity Auios is very small. Thus, we 
suspect that a small line-of-sight velocity goes along with an 
increased probability of misclassihcation (but still less than 
50%). 

To examine this more closely, we plot the fraction of 
correctly classified satellites per line-of-sight velocity bin in 
Figure M The error bars give the standard deviation and 
are quite big for two reasons: hrst, the size of the disadvan¬ 
tageous region of an orbit with given eccentricity (or circu¬ 
larity) can change strongly with the direction of the line of 


sight and secondly, there are only a few satellites with high 
velocities, so there result large deviations in the probabilities 
in the corresponding bins. 

Nevertheless we find the trend sought after: For low 
velocities |Auios| ~ 0, the correct fraction is of the order 
of 60%, whereas for high velocities we nearly gain a 100% 
probability of correct classihcation. An observer could now 
make use of the different probabilities of misclassihcation 
by weighing the observed data with an appropriate func¬ 
tion giving high-velocity data more weight than low-velocity 
data. We applied a simple Auios-weighing to the data used 
for the preparation of FigureEland when averaging the newly 
derived prograde fractions over all nine simulations, we now 
gain a mean of 54.4%. There is still a discrepancy due to 
the non-circularity of the orbits, but the original signal is 
enhanced and closer to the “real” 3D signal. 

A similar dependency can be observed for subhaloes 
with small projected distances from their host halo centre: 
These also are more probable to be misclassihed. We ex¬ 
amined this effect by just neglecting all satellites with a 
projected distance smaller than a specified fraction of the 
virial radius of its host. For higher minimum distances we 
yielded marginally growing observed prograde fractions, up 
to 54.0% when neglecting all satellites within 40% of the 
virial radius of the corresponding host. 

4.5 The effect of interlopers 

Observationally it is quite difficult to determine whether a 
satellite belongs to the host halo or just happens to lie in 
front of (or behind) the host along the line of sight. To model 
these effects and quantify the influence of such interloping 
subhaloes, we included all objects in the 2D analysis that lie 
within the cylindrical tube centred about the position of the 
host halo and with radius Rvii- We only consider satellites 
out to a maximum distance of 5 Rvii (cf. Section 13. 2. 511 as 
possible sources of contamination. 

In Table [7| we summarize the resulting (observed) pro¬ 
grade fraction when including interlopers. These numbers 
are accompanied by the fractional increase of the number of 
satellites and the absolute difference to the situation where 
interlopers are neglected. The new mean fraction is hardly 
affected and only shows a marginal decrease of about 1%. 
We therefore conclude that interlopers have an insignihcant 
effect on the signal as measured by an observer. 

4.6 Alignment of disk and dark matter halo 
angular momentum 

The problem still remains that an observer will only be able 
to measure the distribution and motion of visible matter, 
not the dark matter haloes themselves. Thus, the question 
arises, whether the observed trend of prograde fraction still 
holds for visible matter. The observer determines the rota¬ 
tion axis of a disk galaxy by measuring red- and blueshifts in 
the disk material. We may safely assume that the rotation 
axis of the disk is well aligned with its own angular mo¬ 
mentum, but there are doubts about the alignment of disk 
and dark matter halo angular mo mentum. For example, the 
simple picture of disk formation bv iFall fc EfstathioiJ dlQSfll 
leads to an alignment of disk and halo angular momentum, 
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host 

^int / ^sat 

Tobs 

Jp 


Cl 

1.41 

0.521 

-0.005 

C2 

1.68 

0.521 

-0.004 

C3 

1.40 

0.516 

0.010 

C4 

1.33 

0.601 

-0.008 

C5 

1.32 

0.498 

0.001 

C6 

1.53 

0.531 

-0.015 

C7 

1.45 

0.529 

-0.001 

C8 

1.18 

0.500 

0.001 

G1 

1.41 

0.547 

-0.008 


< fp 0.529 ± 0.029 


Table 7. The interlopers and their effect on prograde fractions, 
riint is the number of satellites including interlopers, A/°^® is 
the absolute difference between the observed prograde fraction 
including interlopers and the previously found value. 

but e.g. lBailin et al. 1^200f^^ found that the disk rotation axis 
and the minor axis of the dark matter halo coincide very well 
within 0.1 Rvii, probably due to their joint evolution. How¬ 
ever, the latter study also emphasizes that the presence of 
baryonic ma tter changes the orientatio n of the dark matter 
halo and in iBailin fc Steinmetd (l200R^ a typical misalign¬ 
ment between the angular momentum vector and the minor 
axis of 25° was found. However, since we have seen in our 
analysis that a change of the host angular momentum only 
had effects of typically a few percent in the prograde frac¬ 
tion, we dare to say that the influence of a disk should not 
change the result in principle. 


5 DISCUSSION 

In this Section we are going to discuss our results and com¬ 
pare them to the findings (and observations) of other groups, 
respectively. 

While we analyzed our 3D data by determining the 
sense of rotation for each subhalo, counting subhaloes on 
prograde orbits and thus determining a prograde fraction 
of 59%, a different method for inv estigating the direc tion of 
rotation of subhaloes was used bv IShaw et alJ (l2005ll . They 
compared the angular momentum of particles in subhaloes 
(rather than subhaloes themselves) to the host’s angular mo¬ 
mentum and found a strong preference of prograde motion 
(cf. Fig. 31 in their study). This method agrees more or 
less with the utilization of cos<l?net, the net orbital angu¬ 
lar momentum of all satellites as presented in Section [3.2.31 
However, applying this method to our data, we also obtain 
a rather large fraction of « 89% corotating (satellite) parti¬ 
cles. 

Our results ar e further consi s tent with the theoreti¬ 
cal signal found by lAzza.ro et a,1. 1 (l2^^5^ , who determined 
the observed prograde fraction for one dark matter halo 
drawn from a cosmological simulation. They used the same 
technique to project the simulation data into the observer’s 
plane, but additionally mapped magnitudes onto their dark 
matter subhaloe s using the (semi-analytical) prescription of 
iKravtsov et a,1. | (l2nn4fl . They found an observed prograde 
fraction of around 55-60%, depending on the mass (or mag¬ 
nitude, correspondingly) of the satellites: massive satellites 


were found to reside preferentially on prograde orbits. We 
could not detect such a mass dependence in our simulations, 
as our mass-weighted curves did not necessarily enhance 
the signal. However, their minimum signal of 55% prograde 
satellites (when including also low-mass satellites, see their 
Tab le 2) agrees we l l with our result of 53%. 

lAzzaro et al. I ^2^^5^ further analyzed a set of obser¬ 
vational data from the SDSS catalogue and actually found 
a strong signal of 61% prograde satellite galaxies. This is 
in contrast to our results: we merely predict a fraction of 
about 53% corotating subhaloes which actually poses a chal¬ 
lenge for an observer to detect. We could imagine several 
reasons for this discrepancy, whose (tiny) individual effects 
ma y add up and thus bri ng our findings into agreement with 
the lAzzaro et al~l ^2005^ results: 

• Determining the mean prog r ade fr action. Using the ob¬ 
servational data, lAzzaro et al. I ll2^^5^ selected 43 similar 
host galaxies and considered the total sample of 76 satellite 
galaxies to belong to a single, fictitious primary. In contrast 
to this, we calculated prograde fractions of subhaloes for ev¬ 
ery individual host halo, before averaging the result. When 
adhering to the same approach of stacking the data prior to 
determining the respective fractions, we end up with 53.9% 
prograde orbits. Grouping together only satellites of host 
haloes which are of roughly the same size and mass (C2, 
C3, C4, C5, C6) leads to the same result. 

• Missing baryonic matter in subhaloes. Most probably 
not all of our subhaloes would carry baryonic matter, some 
of them would therefore escape detection by an observer. It is 
likely that such dark haloes have a low total mass or low cen¬ 
tral density. If we thus restrict our subhaloes and introduce a 
larger minimum mass, minimum density or a lower threshold 
for the maximum value of the rotation curve (Umax, similar 
to observations), we might get a better impression of the 
true satellite galaxy distribution and orientation. We exper¬ 
imented with several threshold values, but were unable to 
obtain mean prograde fractions above 54%. 

• Neglecting satellites with a small projected distance. 
Since it is intrinsically difficult to clearly detect satellite 

S 'es which are very close to the centre. lAzzaro et al. I 
neglected all satellites with a projected distance 
smaller than 20 kpc. This can be advantageous, since satel¬ 
lites close to the centre (in projection) are more likely to 
be misclassified than satellite galaxies residing at farther 
distances (cf. Section 14.4.41 . Thus, the effect of misclassifi- 
cation might decrease, when neglecting such suhaloes with 
small projected distances to the centre of their respective 
host. We found in our simulations only a minor improve¬ 
ment of the signal by at most 1%. 

• Number statist ics. The observational sample of 
lAzzaro et al. I ^2005^ consists of only 76 satellites, whereas 
our total sample contains 1850 subhaloes. Thus, the differ¬ 
ence in the prograde fraction might still be a statistical vari¬ 
ance due to low numbers. 

In summary, we are unable to reproduce such a 
high value of 61% observed p rograde orbits as foun d by 
lAzzaro et al. I ijioO^). However. IZaritskv et al.1 ^1997^ pre¬ 
sented in their observational study 115 satellite galaxies be¬ 
longing to 69 host galaxies (isolated spiral galaxies). Orbit 
orientations were determined for 95 of the satellites with the 
method described in Section r4.2l They found 49 satellites on 
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prograde orbits, yielding a prograde fraction of 51.6%, which 
is perfectly consistent with our prediction of an observed 
pro grade fraction. 

ICarignan et al.1 il99i1l studied the disk galaxy NGC 
5084 with 8 satellites, of which all except one were found 
to be on a retrograde orbit. T his neither supports the pre¬ 
diction of lAzzaro et al. 1 11200, ’ll nor supports our results. An 
excess of retrograde orbits may be explained by the fact that 
dynamical friction is enhanced for satellites orbiting in the 
same direction as the rotation of the host g alaxy. This ef¬ 
fect was studied bv lOuinn fc GoodmarJ (Il986ll using A^-body 
simulations, who found the lif etime of retrograde satellites 
consid erably larger. However, IPenarrubia. Krouna fc Boilvl 
^2002^ found the effect to give only an increase of a few 
per cent in lifetime, depending on the triaxiality of the host 
galaxy halo, so this effect alo ne could not explain the major¬ 
ity of retrograde satellites in ICarignan et ahl Jl997t) . How¬ 
ever, since the number of satellites considered in that study 
is rather low, we consider this signal statistically insignifi¬ 
cant. 


6 CONCLUSION 

We investigated the sense of orientation of satellite orbits in 
cosmological dark matter haloes at redshift z = 0. Our set 
of simulations consisted of eight cluster-sized haloes and one 
simulation of a Milky Way-type galactic dark matter halo. 

From the theory of the generation of angular momen¬ 
tum in dark matter haloes, we expected to get a trend to¬ 
wards satellites corotating with the host. 

We found a mean majority of 59% prograde orbits av¬ 
eraged over all nine simulations. Three of our simulations 
expressed a rather isotropic distribution of orientations of 
orbits with respect to the host. Two of them could be as¬ 
signed a rather small spin parameter of the host, while the 
third could be confirmed as a recent (triple) merger leading 
to random orientations of satellite orbits. 

We validated our results by carrying out a Kolmogorov- 
Smirnov test for the cumulative distribution of angles be¬ 
tween satellite orbit angular momentum and host halo spin. 
All simulations, except the already mentioned exceptions, 
are clearly not in agreement with an isotropic distribution. 
There is also no sign for a symmetry between prograde and 
retrograde orbits in these simulations, stressing that the ex¬ 
cess of corotating satellites is a reliable prediction of the 
hierarchical AGDM structure formation scenario. 

We chose to use the inner host angular momentum, 
i.e. within 30% of the virial mass Mvir (corresponding to 
~ 0.2 7?vir) for defining the rotation of the host halo. How¬ 
ever, other choices (at 50%, 100% of Mvir) lead to the same 
mean prograde fraction of 59%. Thus, uncertainties in the 
angular momentum only have a minor effect. 

We further checked the influence of the environment 
outside the virial radii i?vir of the hosts by including satel¬ 
lites up to 5i?vir. Even though the prograde fraction de¬ 
creased and the orientations became more isotropic, the sig¬ 
nal yet remained - in agreement with the picture that the 
environment of matter and satellites determines the sense of 
host halo rotation. 

Thus, we conclude that there is a trend towards an over¬ 
balance of corotating subhaloes in AGDM simulations at 


redshift z = 0, regardless of the inner host boundary defining 
the host angular momentum and the outer boundary deter¬ 
mining the sphere within which satellites are included. We 
could not confirm any relation between triaxiality and pro¬ 
grade fraction and, even more important, we find our signal 
in the cluster-sized simulations as well as in the galaxy-sized 
simulation. 

The question now remained, how feasible it is to obser- 
vationally test the prediction of having roughly 60% satel¬ 
lites on orbits corotating with their host halo. To obtain an 
answer, we chose 100 random lines of sight projecting our 
three-dimensional data into an observer’s plane. The signal 
still remained, yet weakened to 53% ± 3%, which is difficult 
to be detected. 

The origin of this weakening lies in the observer’s 
method for distinguishing between pro- and retrograde or¬ 
biting satellites by comparing their motion with the part of 
the host where they reside (cf. Section r4.2t . This commonly 
used method works perfectly for circular orbits, yet satellite 
orbits are not necessarily circular but rath er exhibit a distri¬ 
bution of circularities (cf. iGill et al. l2004^ . We showed that 
the recovery of the original signal can be improved by weigh¬ 
ing the satellites with their relative line-of-sight velocities, 
yielding an observed prograde fraction of 54%. 

Additionally, we studied the effect of interlopers, satel¬ 
lites which pass by the observed cluster or galaxy along the 
line of sight, and found it to be negligible (i.e. less than 1%). 

Gomparison wi th observational work y ields different re¬ 
sults: the sturly_c3fJZajitsl g^t al. I (Il997ll agrees with our 
findings, lAzzaro et al. I ll2005^ ’s observati onal results give 
a muc h higher prograde fraction, whereas iGarignan et al. I 
il997ll found for NGG 5084 a strong excess of retrograde 
satellite galaxies. More galaxy-sized simulations, maybe 
even including baryonic matter, would be necessary in order 
to study the prograde fractions on smaller scales in more 
detail, which would then enable a better comparison with 
observational results. 
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APPENDIX A: WEAKENING OF THE 
OBSERVED PROGRADE FRACTION 


Assuming that the fraction of correctly classified satellites 
is the same for satellites on prograde and retrograde orbits, 
we can make following considerations in order to estimate 
the “observed” prograde fraction for our simulations. 

Let Up be the number of prograde and Ur be the num¬ 
ber of retrograde orbits in one of our 3D-simulations. An 
observer of the satellites in this simulation may be mistaken 
in classifying some of the satellites due to problems referred 
to in Section r4.4.1l Thus, the number of prograde orbits 
will be the sum of correctly classified prograde orbits 
and satellites, which are actually retrograde, but have been 
wrongly classified 

_ „right false / . , x 

/ Lp ! bp \ ! bp • y Ai- T j 

Equivalently, the number of retrograde orbits is the sum 
of correctly determined retrograde orbits and wrongly 

classified retrograde orbits 

nf = = -f n“=® (A2) 


Since we know that the latter orbits can only be prograde 
orbits in truth and vice versa, we can write: 


T right 


= rir) — n, 


right 


(A3) 


Let us now assume that the probability of correct clas¬ 
sification c is constant, i.e. independent of the sense of rota¬ 
tion: 


7%p yXf' 


(A4) 


Equation lAlll can then be re-written using equations ifOl . 
llA4ll and we obtain: 


= -b 

= n"*^*'* -b {nr - 

= C - Up + {Ur — C ■ Ur) 

= C ■ Up + (1 — C) ■ Ur 


(A5) 






16 Warnick & Knebe 


For q representing the constant ratio of prograde to retro¬ 
grade orbits in our 3D-simulations, 


q ■.= 



(A6) 


we can express the number of observed prograde orbits by 


Up = c - q • Ur + {I — c) • Ur 
= \c - q +{1 - c)] ■ rir . 


(A7) 


Similarly, we yield for the observed number of retrograde 
orbits: 

= C ■ Hr + {1 — c) ■ Up 

= C ■ Ur + {1 — c) ■ q ■ Ur (AS) 

= \c +q- {1 - c)] ■ rir . 


The ratio of observed prograde to observed retrograde orbits 
is then: 


obs ^ ^ c • g -b (1 - c) 

^ \nr) c + q ■ {1 — c) 


(A9) 


which is not simply the same as q\ 

We conclude that the observed ratio of prograde to 
retrograde orbits is only equal to the real ratio g, if c = 1 or 
c = 0, meaning all or none satellites are classified correctly. 
Note, that the case of 100% prograde (g = 1) or 100% 
retrograde satellites (g = 0) must not be allowed, because 
then our assumption of a constant probability of correct 
classihcation c fails - it will not be well defined anymore. 
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by the author. 



